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This is an overview of pulsar accretion modeling. The physics of pulsar accretion, i.e., the process 
of plasma flow onto the neutron star surface, can be constrained from the spectral properties of 
the X-ray source. We discuss a new implementation of the physical continuum model developed 
by Becker and Wolff (2007, ApJ 654, 435). The model incorporates Comptonized blackbody, 
bremsstrahlung, and cyclotron emission. We discuss preliminary results of applying the new 
tool to the test cases of Suzaku data of CenX-3 and XTEJ1946-I-274. CenX-3 is a persistent 
accreting pulsar with an O-star companion observed during a bright period. XTE J1946-I-274 is 
a transient accreting pulsar with a Be companion observed during a dim period. Both sources 
show spectra that are well described with an empirical Eermi Dirac cutoff power law model. We 
extend the spectral analysis by making the first steps towards a physical description of CenX-3 
and XTEJ1946-f274. 
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Flow of Material comtonized 



Figure 1 Left: Schematic representation of accretion from a wind-emitting star onto a magnetized 
neutron star. After Davidson & Ostriker (1973). Right: Illustration of plasma flowing through the 
accretion column onto the magnetic pole of a neutron star. The soft photons are created through 
blackbody, bremsstrahlung and cyclotron emission. After Becker & Wolff (2007). 

1. Introduction 

X-ray pulsars are binary systems in which material from the donor star is accreted onto a 
neutron star with a strong magnetic field. The first physical continuum model for accreting pulsars 
was provided by two of the coauthors, i.e., the Becker and Wolff 2007 (B&W) model Becker & 
Wolff (2007). The left panel in Figure illustrates how as the material approaches the star, it 
couples with the magnetic field lines at the Alfven radius. As the material falls towards the surface 
of the neutron star, it deccelerates and produces bremsstrahlung emission. Cyclotron emission is 
also present due to the pulsar’s strong magnetic field. Finally, as the material settles on the surface 
of the neutron star, it forms a thermal mound which emits blackbody radiation. The B&W model 
is a bulk and thermal Comptonization model which incorporates the three types of seed photons 
(blackbody, bremsstrahlung, and cyclotron) that are then Comptonized into the approximate power 
law spectrum that we see in accreting pulsars, see right panel in Figure]^. We are working on 
implementing this model for statistical fitting, e.g. in xspec or isis, and making it publicly 
available. Currently phenomenological models are generally used to describe the cutoff power law 
continuum shape of accreting pulsars. 

2. Observations 

We will first introduce the sources, Cen X-3 and XTE J1946-1-274, and discuss empirical mod¬ 
eling of their Suzaku spectra. We will then present the fist results of applying the B&W model to 
these spectra. 

2.1 Centaurus X-3 

Cen X-3 was the first X-ray pulsar discovered in a 1971 observation, with Uhum Giacconi et al. 
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Orbital Phase 
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Figure 2 Left - Suzaku-X\S lightcurve of an observation of a full Cen X-3 orbit in 2008, December. 
The spectral analysis was performed on the data highlighted in red (after Gottlieb et al., in prep.). 
Right - Swift-BAT lightcurve of XTE J1946+274 with the time of the Suzaku observation marked 
as the red vertical line. The upper-left panel shows a closer view of the 2010 series of outbursts 
with apastron and periastron times (after Marcu et al., in prep.). 


(1971). CenX-3 is a persistent High Mass X-ray Binary (HMXB) composed of a neutron star with 
a mass of 1.21 ±0.21 Mq and a pulse period of ~ 4.8 s which orbits an 06.5II-type companion with 
a mass of 20.5 ±O.7M0 Ash et al. (1999). The X-ray source is eclipsed ~ 20% of its 2.1 d orbit. 
Cen X-3 also exhibits several spectral features: three iron (Fe) fluorescent lines at 6.4 keV, 6.7 keV 
and 6.9 keV Naik et al. (2011), a Cyclotron Scattering Resonance Feature (CRSF) at ~ 30keV 
Suchy et al. (2008), and a broad “10 keV bump” at ~ 13 keV, a feature which is relatively common 
in accreting pulsars Coburn et al. (2002); Mihara (1995). 

We studied a 90 ks Suzaku observation spread out over a full orbital period observed in 2008, 
December 8 — 10, eclipse to eclipse, as can be seen in the left panel of Figure^ The spectra shown 
in the left panel of Figure ^ were extracted from a bright part of the lightcurve (shown in red in 
the left panel of Figure ^ in which the hardness ratio was approximately constant. Considering a 
distance to the source of ~ 8kpc Krzeminski (1974), the 1 — 40keV luminosity for the extracted 
bright part is ~ 1.2 X 10^^ erg/s. 

We modeled the XIS and PIN spectra (Gottlieb et al., in prep.) by taking into account the 
following components in xspec: the cross-normalization between the instruments modeled with a 
constant, the partial covering absorption modeled with tbnew_pcf and tbnew_f eo^ Wilms 
et al. (2000), the three Fe emission lines modeled with Gaussian lines (gauss), the 13keV feature 
modeled with another Gaussian line, and a CRSF at ~ 30 keV modeled with a Gaussian optical 
depth profile (gabs)^. The continuum was best fitted with a Fermi Dirac cutoff model Tanaka 
(1986) described by: 


MfdcutiE) E 


1 ±exp 


/ F-Fcut 

V ■S'fold 


-1 


( 2 . 1 ) 


At energy E the photon flux is described by a power law with a photon index F, multiplied by an 


*Tbnew is an updated version of tbabs, see http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/ for more 
details. 

^Three additional Gaussian lines were included to model weak residuals. They will be discussed in Gottlieb et al., 
in prep. 
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Figure 3 Top: Suzaku spectra of the Cen X-3 data highlighted in red in the left panel of Figure 2 
(left) and Suzaku spectra of XTE J1946+274 (right) with the Fermi Dirac cutoff power law best 
fit model (black line). Middle: Residuals of the fitted absorbed continuum without the Fe line 
and CRSF Bottom: Residuals after including the Fe Gaussian line and CRSF Selected best fit 
parameters are listed in Table 1. 


exponential cutoff at energy E^ut with a folding energy Ffoid- The final besf hi model is: 

44cenX-3 (E) = const X tbnew_pcf x tbnew_feo x (power x fdcut 

+3 X gaussFe + gaussiokeV + 3 x gauss^eak) x gabscRSP- (2-2) 

Selecfed parameters of this empirical fit are listed in Table The continuum parameters are high¬ 
lighted. Those are the standard model parameters that describe the continuum which we are trying 
to describe with the physical B&W model. 

Applying an analytical test implementation of the full B&W model (in Fortran), we performed 
a fit “by eye” to the unfolded, unabsorbed Suzaku data of Cen X-3. The result is shown in the left 
panel of Figure^ with the physical parameters listed in Table ^ The results are promising since the 
model qualitatively agrees with the data. The emission is bremsstrahlung dominated overall, with 
a strong contribution from cyclotron emission. The CRSF is not implemented in the model and, 
therefore, we can see the model overpredicting the flux around 30 keV where the observed spectrum 
shows a broad dip due to the magnetic resonant scattering of the photons. The Fe fluorescent lines, 
which are also not part of the continuum model, can be seen in the 6-7 keV region. 

2.2 XTEJ1946+274 

XTE J1946-1-274 was discovered in 1998 September by the All-Sky Monitor (ASM), on board 
of the Rossi X-ray Timing Explorer (RXTE) Smith & Takeshima (1998). In contrast to the bright 
and persistent Cen X-3, XTE J1946-1-274 is a transient HMXB which has shown only two series 
of outbursts: three months in 1998 at discovery and several months in 2010-2011. The optical 
companion is a giant Be-IV-IVe star Verrecchia et al. (2002) orbited by a ~ 15.7 s pulsar Muller 


4 






































































Looking into the Theory of Pulsar Accretion: Cen X-3 and XTE J1946+274 


Diana M. Marcu 


et al. (2012). The orbital period is ~ 169.2 d, the orbital inclination is ~ 46° and the distance to the 
source is 9.5 ±2.9kpc Wilson et al. (2003). The source shows ~ 2 outbursts per orbit, an unusual 
behavior for Be systems in which only one outburst per orbit is generally observed. During the 
1998 series a CRSF was detected at ~ 35 keV Heindl et al. (2001). XTE J1946+274 also shows an 
Fe fluorescent line at 6.4 keV. 

We extracted Suzaku XIS and PIN spectra from a 42 ks observation taken on 2010 October 
12, during a minimum between two outburst from the 2010 series (see right panel of Figure^. 
During the observation the source had a 1 — 40keV luminosity of ~ 3.4 x 10^® erg/s (much lower 
than that of Cen X-3). Similarly to Cen X-3, we modeled the XIS and PIN data (Marcu et al., 
in prep.) taking into consideration the instrument cross calibration, absorption, a Gaussian Fe 
line, and a weak CRSF line. Even though these sources have very different characteristics (e.g. 
luminosity, optical companion, frequency of bright periods, orbital periods, etc.) their continua are 
well described with the same empirical model, fdcut^. The best fit model for XTE J1946-1-274 
in xspec notation is: 

4Txteji 946+274 (-F) = const X tbnew_feo x (power x f dcut-|-gausspe) x gabscRSF- 

(2.3) 

Selected parameters for this fit are listed in Table |I[ The continuum paramteters for XTE J1946-1-274 
are also highlighted. 

Similarly to Cen X-3, we performed a “fit by eye” to the unfolded, unabsorbed Suzaku data of 
XTE J1946-1-274 using the new B&W model implementation. As can be seen in the right panel in 
Figure the continuum is again dominated by bremsstrahlung emission, but with a much weaker 
contribution from cyclotron line emission than we found in Cen X-3. The model again agrees 
well with the data, particularly at lower energies. The physical model parameters applied towards 
this “fit” are listed in Table As expected, due to the difference in luminosity between the two 
sources, the mass accretion rate, electron temperature inside the column, as well as the radius of 
the accretion column are lower for XTE J1946-I-274 than they are for the high-luminosity source 
Cen X-3. 

3. Current and Future Applications 

An earlier application of the B&W model was carried out for the source 4U 0115-1-63 using 
broad-band BeppoSAX data accumulated during an outburst decay, see Ferrigno et al. (2009) and 
xspec model available from C. Ferrigno. In this case the emission from the column is dominated 
by Compton scattered cyclotron emission because of the relatively low magnetic field of this source. 
Its contribution is dominant above ~7 keV. Below this energy, it was necessary to include a thermal 
Comptonization component, which is interpreted as coming from a scattering halo surrounding the 
neutron star surface. 

The “compMag” model Farinelli et al. (2012) provides a numerical solution of the radiative 
transfer equation in a cylindrical column with approximate cross sections for scattering in a mag¬ 
netized plasma using similar assumptions as the B&W model. In the version available for Xspec, 

^Note that this empirical description is not unique, see Marcu et al., in prep, for the application of all commonly 
used empirical models, yielding comparable fit quality. 
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Figure 4 Fit “by eye” applying the B&W physical model to the Suzaku data (black) of Cen X-3 (left) 
and XTE J 1946-1-274 (left). The three types of Comptonized emission are shown: bremsstrahlung 
in green, cyclotron emission in dark purple and back body emission in blue. The physical parame¬ 
ters are listed in Table 3. 

only blackbody seed photons are taken into account in an exponential distribution along the accre¬ 
tion column. This model allows the user to switch between two different velocity profiles of the in 
falling plasma. Work is in progress to include emission from the magnetized optically thin plasma 
of the accretion stream and will be presented in a forthcoming publication, see also Ferrigno et al. 
(2014). 

Here we presented promising first steps towards the physical modeling of accreting pulsars 
using a new analytical implementation of the B&W model. In order to determine the robustness 
of the physical parameters statistical modeling is required and we are currently implementing the 
analytical model into xspec. The analytical and numerical versions of the model will be tested 
on observed X-ray spectra of several accreting pulsars starting with Cen X-3, XTE J1946-I-274, 
4U 0115-1-63 and GX304—1. We plan to conduct a thorough comparison between the empirical 
and physical fits. This will allow us to physically interpret data of accreting pulsars previously 
analyzed with empirical models. 
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CenX-3 XTEJ1946+274 




- 

/ 


- 


1.67(3) 



9 ns+0.04 c 
^■^■^- 0.05 

r 


0.57(2) 

£cut [keV] 


^•^^-0.05 

•Efoid [keV] 

n /\+0.07 
'•^-0.06 

8.9(4) 

£'crsf [keV] 

30.3(9) 

35.2/};^ 

Zrld/dof 

1.31/500 

1.17/466 


Table 1 Main best fit model parameters for the 
Fermi Dirae eutoff power law fits to the Suzaku 
speetra of Cen X-3 and XTE J1946+274. The 
eontinuum parameters are highlighted. The full 
list of fit parameters ean be found in Gottlieb 
et al, in prep, for Cen X-3 and Mareu et al., in 
prep. forXTEJ1946-h274. (a)x lO^^em^^; (b) 
keV^^ em^^s^^; (e)xlO^^. 


Cen X-3 XTE J1946-h274 

M[xl0i’g/s] 

11.32 

1.4 

TeixlO’K] 

3.74 

3.4 

B[x10^2g] 

2.55 

3.04 


700 

100 

^if 

1.01 

2.141 


9.317 

12.79 


Table 2 Physieal parameters used for a “fit by 
eye” of the B&W model for the Suzaku data 
of Cen X-3 and XTE J1946-I-274. The physieal 
parameters are: mass aeeretion rate, M; elee- 
tron temperature inside the eolumn, 7^, mag- 
netie field sfrengfh in fhe emission line region, 
B; radius of fhe aeerefion eolumn, ro; fhe seal¬ 
iering eross seelions parallel lo fhe magnelie 
field (a||) and fhe average eross seefion (a) as 
a funelion of fhe Thompson eross seefion, 07 . 
(d)xlO-^aT. 
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